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4 

Summary 

Boundary  layer  stability,  and  the  ability  to  maintain  a  laminar  flow  near  the  wall,  are 
compromised  in  presence  of  free-stream  vortical  forcing.  The  low-frequency  free-stream 
disturbances  penetrate  boundary  layer  despite  the  sheltering  ability  of  the  mean  shear. 
These  disturbances  form  streaks  via  the  lift-up  mechanism.  The  streaks  subsequently 
become  unstable  and  breakdown  to  turbulence.  The  stability  of  laminar  boundary  layers 
is  examined  in  presence  of  a  thin  wall-film,  of  different  viscosity.  The  influence  of  the 
film  on  the  evolution  of  free-stream  vortical  modes  and  the  generation  of  boundary  layer 
streaks  is  solved  analytically.  It  is  shown  that  a  lower  viscosity  film  can  reduce  the 
transient  amplification  due  to  the  lift-up  mechanism,  and  hence  enhance  the  flow  stability. 
However,  the  presence  of  the  wall-film  also  introduces  a  new  mechanism  for  the  generation 
of  normal  vorticity.  This  mechanism  become  dominant  when  the  viscosity  of  the  lower  film 
is  reduced  below  a  critical  level.  The  results  from  the  linear  analysis  are  complemented  by 
direct  numerical  simulations  which  take  into  account  non-parallel  and  non-linear  effects, 
e.g.  finite  displacement  of  the  two-fluid  interface. 


1.  Introduction 

The  mechanics  of  two-fluid  shear  flows  are  significant  in  many  engineering  and  aero¬ 
dynamical  applications.  For  example,  the  instability  of  two-fluid  flows  affects  the  aero¬ 
dynamic  lift  of  airfoils  in  the  presence  of  deicing  agents,  and  the  heat  transfer  rates  in 
sheared-films  in  cooling  flows.  Here,  the  influence  of  a  different  viscosity  film  near  the 
wall  on  boundary  layer  transition  is  investigated. 

In  single-fluid  boundary  layers,  the  mean  shear  limits  the  penetration  of  vortical  dis¬ 
turbances  towards  the  wall.  This  inviscid  phenomenon  is  known  as  shear- sheltering.  The 
presence  of  viscosity,  however,  renders  the  sheltering  mechanism  less  effective,  and  vor¬ 
tical  disturbances  permeate  the  shear.  As  such,  laminar  boundary  layers  become  prone 
to  bypass  transition  even  at  moderate  levels  of  free-stream  turbulence  (Morkovin  1969). 
In  this  report,  we  focus  on  the  effect  of  introducing  a  thin  wall  film,  with  respect  to  the 
penetration  of  free-stream  vortical  disturbances  towards  the  wall.  The  starting  point  of 
our  study  is  the  continuous  spectrum  eigenmodes  of  the  Orr-Sommerfeld  equation  -  the 
linear  perturbation  equation  for  parallel  shear  flow.  The  dependence  of  mode  shape  on 
wavenumber,  mean  shear,  and  fluid  properties  demonstrate  the  influence  the  wall  film 
on  vortical  penetration  into  the  boundary  layer. 

1.1.  The  continuous  spectrum 

The  first  formal  mathematical  characterization  of  the  continuous  spectrum,  and  the  an¬ 
alytical  expression  of  its  dispersion  relation,  were  presented  by  Grosch  &  Salwen  (1978). 
The  eigenfunctions  of  the  continuous  spectrum  are  oscillatory  in  the  free-stream,  and 
decay  inside  the  boundary  layer  shear,  as  they  approach  the  wall.  An  example  of  Orr- 
Sommerfeld  continuous  modes,  with  different  frequency,  is  shown  in  figure  1.1.  The  vary¬ 
ing  degree  of  penetration,  therefore,  provides  a  link  between  the  free-stream  turbulence 
spectrum  and  the  mean  shear.  Penetration  is  a  function  of  the  perturbation  wavenumber 
(in  all  three  directions)  and  the  flow  properties.  Figure  1.1  shows  the  effect  of  the  tem¬ 
poral  frequency  of  the  mode.  Penetration  in  increased  as  the  modal  frequency  is  reduced 
(while  the  wall- normal  and  spanwise  wavenumbers  are  held  fixed). 

When  a  wall  film  is  introduced,  the  extent  of  modal  penetration  becomes  a  function  of 
the  film  viscosity.  A  sample  eigenfunction  of  the  continuous  modes  for  two-fluid  boundary 
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O-S  mode,  v 

Figure  1.  Effect  of  uj  on  penetration  depth  ( R  =  1000,  kz  —  7r,  ky  =  tt)  , - uj  —  7r, 

-  uj  —  7r/10, - uj  —  7T/100. 


Figure  2.  High  frequency  oscillations  in  mode  shape  for  an  extremely  dense  him. - ,  real  com¬ 
ponent;  - ,  imaginary  component.  kxRer  —  1/100,  ky  —  7r/ 4,  kz  —  7r,  8  b  =  0.3,  =  1000. 

Note,  0  is  plotted  on  separate  scales  in  the  top  and  bottom  fluids 


layers  is  shown  in  figure  2.  The  figure  demonstrates  that  a  low- viscosity  him  changes  the 
wall-normal  wavenumber  of  the  vortical  mode,  inside  the  him  region.  The  scale-conversion 
mechanism  can  enhance  penetration  of  the  vortical  mode  at  low  viscosity.  This  result,  by 
Zaki  &  Saha  (2009)  dehes  the  accepted,  or  common,  wisdom  regarding  shear  sheltering. 

In  order  to  provide  a  quantitative  measure  of  sheltering,  Zaki  &  Saha  (2009)  dehned 
the  penetration  depth, 


where  \(j)\oo  is  the  free-stream  amplitude  of  the  eigenfunction.  This  norm  is  evaluated  over 
a  range  of  him  viscosities,  and  the  results  are  shown  in  figure  3.  The  optimal  viscosity 
ratio,  v*BT  =  [vb/^tT ,  is  clearly  shown.  As  the  viscosity  ratio  is  increased,  the  bottom 
fluid  is  essentially  an  impenetrable  layer  and  the  single-fluid  boundary  layer  results  are 
recovered.  In  the  opposite  limit,  where  the  bottom  fluid  is  less  viscous,  shear  sheltering 
becomes  dominant  and  penetration  is  suppressed.  Therefore,  a  low-viscosity  him  best 
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Figure  3.  Variation  of  penetration  depth  with  viscosity  ratio  where 
kxReT  =  =  0.3,  A  SB  =  0.2 


suited  for  sheltering  the  boundary  layer  from  free-stream  vortical  forcing,  and  delaying 
bypass  transition. 


1.2.  The  amplification  of  streaks 

Beyond  the  initial  stage  of  boundary  layer  receptivity  to  free-stream  turbulence,  elon¬ 
gated  boundary  layer  streaks  are  formed  and  amplify.  Their  growth  has  been  explained 
by  rapid  distortion  theory  (Phillips  1969),  and  the  lift-up  mechanism  (Landahl  1980). 
An  alternative  explanation  is  provided  by  the  Orr-Sommerfeld  and  Squire  equations  for 
linear  perturbations:  A  three-dimensional,  wall-normal  u-disturbance  tilts  the  mean  vor- 
ticity,  and  thus  generates  a  strong  normal  vorticity  disturbance  (Zaki  &  Durbin  2005, 
2006).  At  low  frequencies,  the  generation  of  wall- normal  vorticity  is  equivalent  to  strong 
streamwise  velocity  disturbance,  or  streaks. 

Despite  the  importance  of  the  vorticity  tilting  mechanism  in  the  amplification  of  per¬ 
turbations  in  a  single-fluid  boundary  layer,  it  has  not  been  investigated  in  the  context 
of  two-fluid  boundary  layers.  Furthermore,  the  presence  of  the  interface  in  a  two-fluid 
boundary  layer  can  lead  to  interfacial  instability  mechanisms  (Yih  1967;  Hooper  &  Boyd 
1983;  Hinch  1984;  Hooper  &  Boyd  1987).  In  this  investigation  we  take  advantage  of  the 
thin  layer  effect  which  stabilizes  the  interfacial  mode,  (Renardy  1987;  Charru  &  Hinch 
2000)  when  the  film  viscosity  is  lower  than  the  outer  stream.  The  stabilization  of  the 
interface  mode,  however,  explains  only  the  large  time  asymptotic  behavior  of  the  system. 
In  studies  of  transient  amplification,  even  a  stable  interface  mode  must  be  taken  into 
account  in  evaluating  the  short  time  energy  amplification. 

We  herein  report  on  the  influence  of  wall-films  of  the  amplification  of  streaks  forced  by 
free-stream  vortical  disturbances.  Our  approach  uses  the  linear  perturbation  equations. 
The  solution  to  an  initial  value  problem  which  governs  the  boundary  layer  response  is 
obtained  using  the  eigen-modes  of  the  Orr-Sommerfeld  and  Squire  equations.  Thereafter 
the  boundary  layer  response  to  low-frequency  vortical  forcing  is  evaluated  for  different 
viscosities  of  the  wall  film.  The  findings  based  on  linear  theory  are  complemented  by 
direct  numerical  simulations  (DNS).  The  DNS  takes  into  account  non-parallel  effects 
associated  with  the  downstream  spreading  of  the  boundary  layer.  In  addition,  the  DNS 
also  accounts  for  the  finite  displacement  of  the  two-fluid  interface  -  a  non-linear  effect 
not  captured  by  the  linear  analysis. 

The  background  to  the  influence  of  free-stream  vortical  modes  on  boundary  layers  was 
presented  in  this  section.  The  next  section,  §2,  includes  the  theoretical  formulation  for 
the  base  flow,  the  governing  equations,  and  the  analytical  solution  of  the  initial  value 
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Figure  4.  Schematic  of  a  two-fluid  boundary  layer  under  the  influence  of  a  continuous  mode 

and  an  interface  mode. 

problem.  An  example  evaluation  of  the  solution  of  the  initial  value  problem  is  discussed 
in  §3,  and  comparison  to  DNS  is  given  in  §4  followed  by  the  conclusion  section  (§5). 


2.  Linear  analysis 

2.1.  The  base  flow 

The  two-fluid  boundary  layer  does  not  satisfy  a  similarity  solution.  Therefore,  the  bound¬ 
ary  layer  equations  must  be  marched  downstream  in  order  to  obtain  the  correct  velocity 
profile  at  a  given  streamwise  location.  These  equations,  and  their  solution  procedure, 
were  given  in  Nelson  et  al.  (1995)  and  are  included  below: 

Fd2vF  +  u^F  -  £  (dvF^T  -  =  0,  (2.1) 


where  £  =  j  ,  V  ~  V  ,  and  F  =  (2i/T  J^i/2  with  'h  as  the  streamfunc- 

tion.  The  evolution  of  the  interface  77*  is  described  by  the  kinematic  condition, 


dr] * 


1(V 

Z  \u 


(2.2) 


where  U  and  V  are  the  mean  streamwise  and  wall-normal  velocities  at  the  interface 
respectively.  The  interface  boundary  conditions  are  continuity  of  velocity  and  stress  at 
77  =  77*,  according  to, 

[F\=  0  ;  [dvF}=  0  ;  \^F}=  0, 

where  [.]  denotes  the  change  across  the  interface,  (.)t  —  (. )b •  In  addition,  F(fl,r])  must 
satisfy  the  boundary  conditions, 


F& 0)  =  0,  dr,F&  0)  =  0,  lim  dr,F&  77)  =  1. 

?7— >00 

The  solution  procedure  of  the  above  equations  follows  the  description  in  Schlichting 
(1987,  pp  187-191),  and  was  implemented  by  Zaki  &  Saha  (2009).  This  solution  provides 
the  mean  profile  used  in  the  computation  of  the  continuous  spectrum  modes. 
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2.2.  The  linear  perturbation  equations 

The  evolution  of  linear  perturbations  in  a  parallel  shear  flow  can  be  fully  described  by 
three  equations  for  the  normal  velocity  u,  the  interface  deflection  /,  and  normal-vorticity 
rj.  These  equations  are  expressed  in  operator  form  according  to, 


d_ 

dt 


0  0  \ 
0 

0  yj 


Under  the  parallel-shear-flow  approximation,  the  equations  are  homogeneous  in  time,  as 
well  as  the  streamwise  x  and  spanwise  2  directions.  Therefore,  the  following  form  of  the 
solution  can  be  assumed, 


(h(v) 

\xj(y) 


(2.3) 


Substituting  the  normal  mode  assumption  2.3  in  the  governing  equation  2.3  yields  the 
Orr-Sommerfeld  ,  interface,  and  Squire  eigenvalue  problems  for  two-fluid  flows, 


-UJ 


/  d2  —  k2  0  0^ 
0  1  0 


0 


0 


'fa(y)\ 

( 

0 

°\ 

(hiv) 

f  = 

-  r «(»-»/) 

0 

f 

.Xj(y )/ 

v 

0 

\Xj(y) 

(2.4) 


The  operators,  j£fj,  ^ ,  ST j  are  defined  as, 

Xj  =  ikxdyUj  +  -T_(d2  -  fc2)2  _  ikxUj(d2  _ 

pOO 

&(•)  =  /  (•)%  -  Vf)dy 

JO 

y  =  -ikxU{yf) 

yJ  =  1T(dl-k2)-ikxuj 

=  ikzdyUj , 


where  k2  =  k2  +  k2  and  dy  =  . 

Boundary  conditions  must  be  enforced  at  the  wall,  across  the  interface,  and  in  the 
free-stream.  The  wall  and  free-stream  boundary  conditions  are  similar  to  the  single-fluid 
case, 


0s(O)  =  O;  dy4>B(  0)  =  0;  Xb(0)  =  0 

<h(y  ^  00)  =  0;  dy4)T{y  oo)  =  0;  xt(v  ^  00)  =  0. 


Across  the  interface,  the  following  conditions  of  velocity  and  stress  continuity  must  be 
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enforced  on  (j)  at  y  =  yf, 


Pt 


=  Pb 


<t>T  =  4>B 

dy(pT  —  ikxdyUrf  =  dy(j)B  —  ikxdyUsf 
!*r(dy  +  k2)4>T  ~  ikxnTd2yUTf  =  HB(d2y  +  k2)4>B  ~  ikxHBd2yUBf 

—  iu)dy<j)T  —  R&a  i  ( dy  —  2>k2\(j)T  +  ikx{UTdy(j>T  —  4> TdyUx ) 


3  \  v 
-i  ( 


-  iujdy(j)B  -  Re  a  I  dy  -  3k  )<f>B  +  ikx(UBdy4>B  -  (j)BdyUB) 


+  We~1kAf. 


The  interfacial  conditions  on  x  are  as  follows, 

XT  ~XB  =  —ikz(dyUT  ~  dyUB)f 
UrdyXT  —  I^BdyXs  =  ikz(iiTd2UT  —  HBd2UB)f 


2.3.  The  solution  of  the  initial  value  problem 

Here,  we  consider  the  effect  of  a  vertical  velocity,  v,  perturbation  on  the  two-fluid  bound¬ 
ary  layer.  We  prescribe  the  initial  condition  for  v,  and  to  compute  the  boundary  layer 
response  using  the  perturbation  equations  2.3. 

For  a  general  disturbance,  we  introduce  the  Fourier  transform  in  the  streamwise  and 
the  spanwise  directions, 

/OO  BOO 

/  Vj(x,t)e~lkxXe~lkzZdxdz 

-oo  J  —  OO 

/OO  BOO 

/  f(x,z,t)e-ik*xe-ik*zdxdz 

-OO  J  —  oo 

/OO  BOO 

/  r]j(x,t)e~lkxXe~lkxZdxdz. 

-oo  J  —  oo 

The  normal  velocity  and  the  interface  equations  form  an  autonomous  sub-system, 
independent  of  the  normal  vorticity.  Therefore,  this  subsystem  is  solved  first.  The  solution 
is  expressed  in  terms  of  the  discrete  and  continuous  eigenfunctions  of  the  Orr-Sommerfeld 
and  interface  equations.  The  discrete  and  continuous  spectra,  combined,  form  a  complete 
basis.  Therefore,  the  solution  is  expanded  according  to, 


The  orthogonality  of  the  adjoint  system  is  invoked  in  order  to  solve  for  the  evolution  of 
the  coefficients,  a £  and  a £  k  .  Here,  we  focus  on  an  initial  disturbance  that  is  composed 
of  only  a  single  Orr-Sommerfeld  mode, 

a£,m(0)  =  0  a£,/cy(0)  =  AS(ky  —  ky ). 

Following  the  procedure  in  Salwen  &  Grosch  (1981),  the  following  solution  is  obtained, 

fk(t)  =  Ah,kye~lWn''‘vt  • 


(2.6) 
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The  Squire  equation 

Zaki  &  Durbin  (2005)  proposed  the  solution  to  the  forced  Squire  equation  for  a  single¬ 
fluid  boundary  layer,  in  terms  of  the  eigenfunctions  of  the  Squire  modes.  The  Squire 
eigenfunctions,  %,  satisfy  the  homogeneous  Squire  equation, 

Hw  +  ikxUj)(d2y  -  k2)xj  ~  -T-(dy  -  k2)xj  =  0  (2.7) 

and  homogeneous  interfacial  conditions, 

Xt  ~Xb  =  0 

TrdyXr  ~  l^sdyXB  =  0.  (2.8) 

For  a  two-fluid  boundary  layer,  the  homogeneous  Squire  eigenfunctions  remain  contin¬ 
uous  across  the  interface.  This  introduces  a  difficulty  since  the  normal  vorticity  we  wish 
to  represent  is  discontinuous  at  the  interface,  and  must  satisfy  the  following  inter  facial 
jump  conditions, 


-^B,k(yf^)  =  -ikz[dyUT{yf)  -  dyUsiVf) J  f%(t), 

fJ,Tdyr]T^(yf,t)  -  yBdyr]Bj:(yf,t )  =  -ikz(^iTd2yUT{yf)  ~  gBd2yUB{yf)^jf^{t). 

In  order  to  satisfy  the  above  jump  conditions  implicitly,  a  change  of  variables  is  intro¬ 
duced.  The  new  variable  T  is  continuous  and  defined  as, 


t)  =  Vtfiv,  t)  +  ikzdyUjfR(t)g(y).  (2.9) 

The  function  g(y)  must  be  continuous  and  differentiable  over  the  interval  [0,  oo)  for  T 
to  be  continuous  and  satisfy  the  interfacial  jump  conditions.  Moreover  g(y)  should  equal 
unity  at  the  interface  and  zero  at  the  wall.  The  choice  of  g(y)  does  not  affect  our  results 
and  any  g(y)  which  satisfies  the  above  conditions  can  serve  the  purpose  of  solving  the 
initial  value  problem.  In  this  particular  case,  we  use  g(y)  =  -^eyf~y .  Thus  the  interfacial 
boundary  conditions  on  T  are  identical  to  those  on  the  eigenfunctions  y, 


Hrdy^T^{yf,t)  -  fj,Bdy^B^(yf,t)  =  0. 

Substituting  definition  2.9  into  the  Squire  equation  yields  the  evolution  equation  for  T, 

JU  .  g(2/,  t)  -  ./,'I'.;:(//.  t)  =  (2.10) 

where  the  forcing  term  J ^  is, 


TPd  -  —ikzA 


dyUj 


huM  ~  huM)9(y) )  +  Rej 1 


^ j,k,k. 


dvUjg(y) 


Following  Zaki  &  Durbin  (2005),  we  expand  T  in  terms  of  the  homogeneous  Squire 
eigenfunctions, 

Ns(lk  n 

h^Xk^dky.  (2.11) 

m=  1  Jky 

The  amplitude  functions,  m(t)  and  k  (t)  can  be  derived  using  the  bi-orthogonality 
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condition  satisfied  by  the  adjoint  Squire  eigenfunctions.  Following  a  similar  procedure  as 
outline  in  Zaki  &  Durbin  (2005)  yields  the  evolution  equations  for  the  spectral  coefficients, 


^Uk,nbk,n^  + 
~tljJk,kybk,ky^)  + 


k  ,  ky  ^ 

iUJk,kyl 


(2.12) 


Equation  2.12  is  a  first  order,  ordinary  differential  equation  and  can  be  solved  exactly 
for  the  amplitude  functions, 


Since  the  continuous  spectrum  of  the  Squire  and  the  Orr-Sommerfeld  equation  overlap, 
there  exists  a  continuous  Squire  mode  which  resonates  with  the  forcing  Orr-Sommerfeld 
mode.  Therefore,  for  the  resonant  mode,  ky  =  ky,  the  amplitude  is  given  by 


k  k  ^  — 2CJ  r  t 

p  /v ,  —  p  k ,  m 


lLJk,ky  liJk,m 


h.kM  =  k.kS°)e  lulR,kyt  ~ 


'k,ky\°/  —  k,k 

and  for  non-resonant  Squire  modes  ky  ky 


(2.13) 


k.kA1)  =  k.kA°)e  ™%'kyt  -  <xc , 


Ky  _  g 


^^k  k  ^^k  k 

r\j  ^  r\j  y  tv  j  i\j  y 


"k,ky\~y  ~k,kyy~y~  ^k,ky 

Finally,  the  normal  vorticity  r]  is  recovered  from  T  using  equation  2.9. 


3.  Evaluation  of  the  boundary  layer  response  to  vortical  forcing 

Although  the  solution  to  the  initial  value  problem  is  general,  our  objective  is  to  evalu¬ 
ate  the  effect  of  viscosity  stratification  on  a  two-fluid  boundary  layer.  Therefore,  surface 
tension  effects  are  ignored,  and  we  assume  unit  density  ratio.  Focus  is  placed  on  study¬ 
ing  the  two-fluid  boundary  layer  response  to  an  initial  forcing  by  a  streamwise  oriented 
vortex.  This  type  of  forcing  is  known  to  be  most  effective  at  generating  streaks  in  the 
single-fluid  problem.  Furthermore,  the  initial  condition  assumes  zero  interface  displace¬ 
ment.  This  condition  is  constructed  from  a  superposition  of  an  O-S  mode  with  finite 
interfacial  displacement  and  an  interface  mode. 

The  continuous  Orr-Sommerfeld  mode  resonantly  forces  Squire  equation,  and  the  so¬ 
lution  2.13  explicitly  demonstrates  initial  algebraic  growth  with  time  for  the  resonant 
Squire  eigenmode.  The  inner  product,  |(xjl  ,  is  shown  in  figure  5,  for  forcing  by 

(i)  the  continuous  Orr-Sommerfeld  mode  and  (ii)  the  interface  mode.  The  peak  in  the 
projection  is  at  the  resonant  Squire  mode,  and  hence  algebraic  transient  growth  is  ex¬ 
pected  in  the  response. 

The  normal  velocity  profile  is  shown  in  figure  6(a)  at  various  times.  It  is  oscillatory 
in  the  free-stream  and  decays  inside  the  boundary  layer  similar  to  the  single-fluid  case. 
The  entire  normal  velocity  field  decays  with  time  owing  to  viscous  dissipation.  At  large 
times  there  is  no  vertical  velocity  present  because  the  continuous  Orr-Sommerfeld  mode 
is  exponentially  stable,  and  there  is  no  vertical  velocity  perturbation  associated  with  the 
interface  mode. 

The  time  evolution  of  the  interface  displacement  is  shown  in  figure  6(b),  and  can  be 
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Figure  5.  Projection  of  forcing  \k\  =  \{x\  k  ->&)\  on  the  homogeneous  Squire  eigenfunc¬ 
tions.  - ,  interfacial  mode; - ,  continuous  Orr-Sommerfeld  mode.  The  forcing  mode  has 

Mint  =  0.000373,  kx  =  0.001,  ky,Cont  =  7r,  kz  =  n  and  the  Reynolds  number  is  Re  —  800.  The 
film  thickness  and  viscosity  are  d  —  0.1  and  hbt  =0.3. 


Figure  6.  (a)  Normal  velocity  and  (b)  interface  displacement  at  different  instants  of  time. 

- ,  t  —  0; - ,  t  =  44; - ,  t  —  314.  The  forcing  mode  has  Ccw  —  0.000373,  kx  —  0.001, 

ky,cont  —  7T,  kz  —  7r  and  the  Reynolds  number  is  Re  =  800.  The  film  thickness  and  viscosity  are 
d  —  0.1  and  hbt  —  0.3. 


explained  by  examining  equation  3.1, 


/&W  —  Ak~kyfk,ky 


-iujt  r  t 


Ak.in.t  f 1 


k,intJ  k,int 


(3.1) 


The  two  amplitudes  and  .  in  equation  3.1  were  chosen  such  that  the  initial 

interface  displacement  is  zero.  However  this  cancellation  is  no  longer  maintained  at  larger 
times  since  the  continuous  Orr-Sommerfeld  mode  decays  much  faster  than  the  interfacial 
mode.  The  net  interface  displacement  therefore  increases  in  time.  At  large  times,  only  the 
interface  mode  is  present  and  therefore  interfacial  displacement  should  decay  as  . 

Figure  7(a)  shows  the  normal  vorticity  response  as  a  result  of  the  forcing  by  the 
interfacial  and  continuous  Orr-Sommerfeld  mode.  The  free-stream  solution  is  trivial  at 
all  instants  of  time.  Normal  vorticity  is  generated  by  the  vorticity  tilting  mechanism, 
whereby  the  normal  velocity  component  of  the  perturbation  tilts  the  mean  shear  and 
generates  perturbation  vorticity.  Since  there  is  no  mean  shear  in  the  free-stream,  no 
normal  vorticity  is  generated.  However  inside  the  boundary  layer,  the  normal  vorticity 
perturbation  amplifies  due  to  the  vorticity  tilting  mechanism. 

Viscous  decay  sets  in  on  a  time-scale,  tviscous  ~  0{Re/{k \  +  ky  +  k%)).  At  larger 
times,  the  normal  velocity  perturbation  vanishes,  and  the  normal  vorticity  equation  is 
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Figure  7.  (a)  Normal  vorticity  response  at  different  instants  in  time, - ,  t  —  0; - ,  t  =  44; 

—  -  — ,  t  —  314.  A  zoomed  in  view  of  the  boundary  layer  is  shown  in  pane  (b).  The  forcing 
mode  parameters  are  ujint  —  0.000373,  kx  —  0.001,  kyjCont  =  7r,  kz  =  n  and  the  Reynolds  number 
is  Re  =  800.  The  film  thickness  and  viscosity  are  d  —  0.1  and  /hbt  —  0.3. 


Figure  8.  Squire  response  in  a  two- fluid  boundary  layer  at  the  time  of  maximum  amplification 

in  energy,  d  —  0.1  - ,  (jlbt  —  1.0; - ,  hbt  —  0.5  ; - ,  hbt  —  0.3; - ,  hbt  —  0.2. 

kx  -  0.001,  ky  =  7T,  kz  —  7T,  Re  —  800 

driven  purely  by  the  displacement  of  the  interface.  In  the  long-time  limit,  a  jump  in 
normal  vorticity  is  observed  across  the  interface  and  the  vorticity  perturbation  decays 
away  from  the  interface,  as  shown  in  figure  7(b).  The  large-time  normal  vorticity  field  is 
identical  to  the  particular  vorticity  associated  with  the  interface  mode. 

3.1.  Viscosity  stratification  effects 

A  reduction  in  viscosity  ratio,  /ibt,  increases  the  amplitude  of  the  interface  mode  and 
hence  the  neutral  interfacial  wave  generates  large  amplitude  normal  vorticity.  The  effect 
of  viscosity  stratification  on  the  vorticity  tilting  mechanism  must  also  be  evaluated. 
Figure  8  shows  the  Squire  response  at  the  time  of  maximum  amplification.  The  peak  in 
normal  vorticity  response  associated  with  the  vorticity  tilting  mechanism  reduces  as  the 
viscosity  ratio  /abt  is  reduced.  The  forcing  term,  which  is  shown  in  figure  9(a-b)  is  also 
attenuated.  Therefore,  it  can  be  concluded  that  the  temporal  amplification  due  to  the 
vorticity  tilting  mechanism  is  reduced  by  lowering  of  Hbt- 

The  origin  of  the  two  peaks  in  normal  vorticity  has  been  addressed;  the  first  maximum 
is  due  to  the  tilting  of  mean  vorticity  and  the  second  is  associated  with  the  neutral 
interfacial  wave.  The  first  peak  is  therefore  dependent  on  the  strength  of  the  forcing 
term,  ^  (Zaki  &  Durbin  2005).  This  term  is  affected  by  the  mean  shear  distribution 
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and  the  variation  in  </>  with  viscosity  stratification.  As  hbt  is  reduced,  the  shear  in  the 
bottom  fluid  increases  and  the  shear  in  the  top  fluid  is  reduced  in  order  to  maintain  the 
mean  shear  stress  continuity  at  the  interface.  The  shape  of  0  in  a  two-fluid  boundary 
layer  was  discussed  in  detail  in  the  work  of  Zaki  &  Saha  (2009).  It  was  reported  that 
the  penetration  of  normal  velocity  perturbation  into  the  boundary  layer  is  inhibited  for 
low  viscosity  ratios  due  to  the  presence  of  stronger  shear  in  the  bottom  fluid.  Therefore 
the  overall  forcing  term  ^  is  reduced  in  strength  as  the  viscosity  ratio  is  reduced.  The 

coupling  coefficient  6  =  ^ proposed  by  Zaki  &  Durbin  (2005)  was  computed, 

and  is  shown  in  figures  10(a-c).  As  the  viscosity  ratio  is  reduced,  the  optimum  coupling 
coefficient  is  significantly  reduced  compared  to  the  single-fluid  value. 

Below  a  critical  viscosity  ratio,  the  neutral  inter  facial  wave  generates  normal  vorticity 
which  exceeds  the  mean- vorticity  tilting  mechanism.  For  a  given  viscosity  ratio  a  thicker 
film  absorbs  the  mean  shear  to  a  greater  extent  and  therefore  the  coupling  between  mean 
shear  and  normal  velocity  is  weaker.  This  implies  that  the  vorticity  tilting  mechanism 
is  expected  to  be  weaker  for  a  thicker  film.  This  prediction  is  substantiated  by  figure  11 
which  shows  the  coupling  coefficient  optimized  over  all  kXyky  and  kz.  The  optimal  0  is 
larger  for  the  thinner  film. 

The  solution  to  the  initial  value  problem  indicates  that  the  temporal  amplification  of 
streaks  is  reduced  in  presence  of  low-viscisoty  wall  film.  The  equivalent  spatial  problem 
is  now  studied  using  direct  numerical  simulations.  In  the  next  section,  the  numerical 
method  used  in  our  direct  numerical  simulations  is  introduced,  followed  by  a  discussion 
of  the  results  from  our  direct  computations. 


4.  Direct  numerical  simulations 

Direct  numerical  simulations  (DNS)  take  into  account  the  downstream  spreading  of 
the  base  flow,  and  non-linear  effects  such  as  the  finite  displacement  of  the  interface  be¬ 
tween  the  two  fluids.  These  phenomena  were  not  modelled  in  the  linear  analyses  above. 
Non-parallel  effects  can  play  an  important  role  because  the  free-stream  vortical  forcing 
of  interest  is  low-frequency.  Therefore,  the  mean  flow  spreading  with  distance  from  the 
leading  edge  effects  the  boundary  layer  response  to  these  elongated  modes.  With  re¬ 
gards  to  non-linear  effects,  finite  displacement  of  the  interface  can  change  the  effective 
amplification  rate  of  disturbances  Cheung  &  Zaki  (2010),  and  ultimately  change  in  the 
interface  topology  Cheung  &  Zaki  (2011).  These  effects  have  modelled  herein  using  direct 
solutions  of  the  Navier-Stokes  equations,  and  the  two-fluid  interface  is  simulated  using  a 
mass-conserving  interface-capturing  scheme. 

In  order  to  compare  the  DNS  results  to  linear  theory,  the  numerical  simulations  are  per¬ 
formed  at  small- amplitude  inflow  perturbations.  A  single  low-frequency  continuous  Orr- 
Sommerfeld  mode  is  prescribed  at  the  inflow  of  our  computational  domain,  and  its  spatial 
evolution  is  computed  for  six  different  viscosity  ratios,  (Ibt  =  {0.2,  0.3,  0.5,  0.7,  0.8, 1.0}. 
The  film-thickness  is  set  to  10%  of  the  boundary  layer  thickness  at  the  inlet. 

4.1.  The  numerical  method 

The  time-dependent  Navier-Stokes  equations  are  solved  using  a  fractional  step  algorithm. 
A  local  volume-flux  formulation  is  utilized,  on  a  staggered  computational  grid  in  general 
curvilinear  coordinates  (Rosenfeld  et  al  1991).  The  convective  terms  are  treated  explic¬ 
itly  using  the  Adams-Bashforth  scheme.  The  pressure  and  diffusion  terms  are  treated 
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Figure  12.  Vertical  profiles  of  streamwise  mean  velocities  at  various  x-stations. 
Hbt  =  1.0; - ,  hbt  —  0.5; - ,  iibt  =  0.3; - ,  hbt  —  0.2. 


by  implicit  Euler  method,  and  by  Crank-Nicolson  method,  respectively.  This  code  is 
parallelized  using  Message  Passing  Interface  (MPI). 

The  inflow  to  the  computational  domain  is  a  superposition  of  the  two-fluid  bound¬ 
ary  layer  profile,  and  a  single  mode  perturbation.  The  inlet  Reynolds  number  is  Res0  = 
800,  based  on  the  inlet  boundary-layer  thickness  and  free-stream  velocity.  The  stream- 
wise,  wall-normal  and  spanwise  sizes  of  the  computational  domain  are,  Lx/5q  =  300, 
Ly/So  =  20  and  Lz/S o  =  8,  respectively.  The  number  of  grid  points  used  is  (Nx,  Ny , 
Nz)=(257, 129,  65).  The  grid  is  uniformly  spaced  in  the  streamwise  x  and  spanwise  2 
directions,  and  is  clustered  near  the  interface  in  the  wall-normlal  y  direction.  The  height 
of  the  computational  domain  is  selected  to  be  sufficiently  large  in  order  to  ensure  that  the 
boundary  layer  is  unaffected  by  the  far  field  boundary  conditions.  The  convective  outflow 
condition  is  applied  at  the  exit  plain.  A  no-slip  boundary  condition  is  imposed  at  the 
solid  wall.  In  the  free  stream,  the  following  boundary  conditions  are  enforced,  u  =  Uoo 
and  —  0.  Periodicity  is  enforced  in  the  spanwise  direction.  For  the  single- mode 

perturbation,  the  inflow  mode  is  uj  =  0.01,  ky  =  7r,  kz  =  7 r,  and  the  amplitude  of  the 
inflow  disturbance  is  normalized  such  that  vrms  =  0.01%  in  the  free  stream. 

The  two-fluid  interface  is  represented  by  an  iso-surface,  or  a  level-set  of  a  smooth 
function.  The  conventional  level-set  technique  uses  the  signed  distance  function.  However, 
this  approach  yields  poor  mass  conservation.  Therefore,  the  current  work  adopts  the 
conservative  level-set  method  in  conjunction  with  a  ghost  fluid  approach  (Desjardins 
et  al  2008).  This  method  ensures  accurate  and  robust  interface  transport  by  using  a 
hyperbolic  tangent  level  set  function.  Moreover,  the  sharp  discontinuity  of  the  interface 
can  be  captured  by  the  ghost  fluid  method. 

4.2.  Numerical  results 

The  influence  of  the  viscosity  ratio  on  the  mean  velocities  is  shown  in  figure  12.  Here, 
the  y- axis  is  normalized  by  the  local  boundary  layer  thickness  and  four  different  viscosity 
ratios  are  plotted.  It  is  clearly  seen  that  the  wall  film  absorbs  the  mean  shear  and  this 
tendency  is  more  pronounced  with  decreasing  viscosity  of  the  film.  As  a  result,  it  is 
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Figure  13.  Vertical  profiles  of  urms  at  various  x-stations. - ,  /ibt  —  1.0; - ,  /ibt  —  0.5; 

- ,  /ibt  —  0.3; - ,  ii bt  —  0.2. 


expected  that  the  lift-up  mechanism  will  be  weaker  for  a  less  viscous  film,  since  the  mean 
mean  shear  is  decreased  in  the  top  fluid  as  shown  in  figure  9. 

The  root-mean-square  (rms)  streamwise  velocity  fluctuations  at  different  ^-locations 
are  displayed  in  figure  13.  The  wall-normal  axis  is  shown  in  the  range  0.3  <  y/Sgg  <  0.7, 
outside  the  wall  film,  in  order  to  demonstrate  the  attenuation  of  the  streak  clearly.  The 
streamwise  velocity  fluctuations  for  /ibt  =  1-0  are  increasing  from  the  inlet  to  x/6g  ~  45 
and  then  decay  as  they  convect  further  downstream.  This  trend  is  also  observed  for  other 
viscosity  ratios.  The  boundary  layer  streaks  are  most  amplified  at  x/5g  «  45,  which 
is  very  close  to  the  time  of  maximum  amplification  predicted  by  the  linear  analysis. 
Furthermore,  the  streak  amplitude  becomes  weaker  with  decreasing  viscosity  ratio.  This 
is  consistent  with  the  linear  prediction  that  lower  viscosity  wall  films  attenuate  the  energy 
amplification  inside  the  boundary  layer. 

The  maximum  urms  is  extracted  at  every  downstream  location  and  is  plotted  in  figure 
14.  In  figure  14(a),  the  maximum  is  computed  from  the  disturbance  field  in  the  outer 
fluid,  and  therefore  reflects  the  strength  of  the  streaks  in  the  outer  stream  due  to  tilting  of 
mean  vorticity.  The  maximum  values  of  the  urrns  for  the  streaks  are  located  at  x/5g  ^40. 
This  peak  amplification  decreases  at  lower  viscosity  ratios,  which  is  consistent  with  the 
prediction  of  linear  theory.  The  relative  decrease  in  the  maximum  urrns  is  as  large  as  30% 
for  the  lowest  viscosity  ratio,  /ibt  =  0.2.  It  should  be  noted  in  figure  14(b),  however, 
that  the  maximum  values  of  urms  for  /ibt  =  0.2  and  0.3  are  significantly  higher  than 
those  of  other  viscosity  ratios  and  the  maxima  show  slow  decay  rates  for  x/Sq  >  100. 
This  indicates  that  the  inflow  disturbances  decay  much  faster  than  the  R-perturbation 
induced  by  the  interfacial  mechanism,  as  explained  by  linear  theory. 

Plan  views  of  the  instantaneous  streamwise  velocity  fluctuations  are  shown  in  figure  15. 
The  figure  also  affirms  the  reduction  of  the  streak  strength  by  a  lower-viscosity  wall-film. 
The  solution  to  the  initial  value  problem  presented  earlier  neglects  both  non-parallel 
and  nonlinear  effects.  Here,  the  base-ffow  spreading  and  the  interfacial  displacement 
are  fully  represented,  but  the  amplitude  of  the  disturbance  was  kept  small  in  order  to 
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Figure  14.  Maximum  urms  versus  downstream  distance  from  the  inlet  plane  for  (a)  the  outer 

fluid  and  (b)  for  the  entire  y-domain. -  hbt  —  1; - I^bt  *=  0.8; - hbt  —  0.7; - 

Iibt  —  0.5;  — B —  /ibt  —  0.3;  — V —  hbt  —  0.2. 


Figure  15.  Instantaneous  contours  of  streamwise  velocity  fluctuations  evaluated  at  y/6 o  =  0.7. 
(a)  iibt  —  0.7,  (b)  hbt  =  0.5,  (c)  iibt  —  0.3.  Contour  levels  correspond  to  —0.002  <u<  0.002 
at  an  increment  of  0.0004. 


preserve  linear  behavior.  The  agreement  between  the  DNS  and  linear  theory  therefore 
complements  the  analyses  presented  herein. 


5.  Conclusions 

The  amplification  of  streaks,  due  to  low-frequency  vortical  disturbances  has  been  in¬ 
vestigated  in  two-fluid  boundary  layers.  The  presence  of  a  different  viscosity  wall- film 
influences  (i)  the  penetration  of  vortical  disturbances  into  the  mean  shear  and  (ii)  and 
the  amplification  of  the  boundary  layer  response.  The  latter  was  studied  by  solving  the 
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initial  value  problem  of  normal- vorticity  generation  (Squire  response)  by  vertical- velocity 
forcing  (O-S  mode).  The  analytical  solution  shows  that  the  presence  of  a  low- viscosity 
wall-film  can  weaken  the  lift-up  mechanism.  This  is  due  to  the  redistribution  of  the  mean 
shear,  which  becomes  concentrated  near  the  wall,  where  the  vertical  velocity  is  vanish¬ 
ingly  small.  This  explanation  was  supported  by  evaluation  of  the  coupling  coefficient 
between  the  O-S  forcing  and  the  Squire  response.  The  coefficient  affirmed  that  the  re¬ 
duction  in  viscosity  ratio  yields  weaker  streaks  in  the  context  of  the  lift-up  mechanism. 
However,  the  presence  of  a  viscosity  mismatch  at  the  interface  introduces  another  source 
of  normal  vorticity.  The  interface  deformation  also  forces  the  Squire  equation,  and  can 
lead  to  normal-vorticity  response,  although  at  a  longer  time  scale. 

The  linear  analyses,  while  shed  light  on  the  linear  mechanism  of  streak  amplification, 
do  not  address  (a)  non-paralllelism  in  the  base  flow  and  (b)  non-linear  effects  such  as 
the  finite  deformation  of  the  interface.  The  former  is  important  because  the  streaks 
are  streamwise  elongated  structures,  and  can  therefore  be  sensitive  to  the  downstream 
spreading  of  the  boundary  layer.  Since  the  interface  displacement  is  finite,  and  Klebanoff 
streaks  reach  high  amplitude,  non-linear  effects  can  be  appreciable.  In  order  to  address 
non-parallelism  of  the  base  flow  and  non-linear  effects  due  to  the  interface  movement, 
direct  numerical  simulations  were  carried  out,  albeit  for  low- amplitude  disturbances.  An 
accurate,  conservative  level-set  approach  was  used  in  order  to  track  the  evolution  of  the 
interface,  and  the  boundary  layer  response  to  particular  vortical  modes.  The  results  were 
consistent  with  the  prediction  of  linear  theory:  the  outer  streaks  were  weaker  in  the 
presence  of  low-viscosity  wall  film.  In  addition,  a  near-interface  peak  in  ^-perturbation 
was  observed  due  to  the  interface  deformation  mechanism.  This  peak  increased  with  lower 
film  viscosity,  and  acted  on  a  longer  time-scale  than  the  well-known  lift-up  mechanism. 
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